Monte Carlo simulation of the electron transport through thin slabs: 
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The Monte Carlo simulation of the electron transport through thin slabs is studied with five 
f"*^ | general purpose codes: penelope, GEANT3, GEANT4, EGSnrc and mcnpx. The different material 

foils analyzed in the old experiments of Kulchitsky and Latyshev [Phys. Rev. 61 (1942) 254-266] 
and Hanson et al. [Phys. Rev. 84 (1951) 634-637] are used to perform the comparison between 
the Monte Carlo codes. Non-negligible differences are observed in the angular distributions of the 
, transmitted electrons obtained with the some of the codes. The experimental data are reasonably 

O |i well described by EGSnrc, penelope (v. 2005) and GEANT4. A general good agreement is found for 

EGSnrc and GEANT4 in all the cases analyzed. 



I. INTRODUCTION 



In Monte Carlo (MC) simulation of electron transport, the simulation of all the interactions suffered by the electrons 
when they go through a material, the so-called detailed simulation (DS), is, in general, unpractical because of the 
long computing times required to reduce the electron energies below the absorption threshold. Therefore, DS is done 
in practice only if the energy of the electrons is low and/or the targets are thin. 

To solve the situation, "condensed history" schemes are employed in much of the usual MC radiation transport 
codes of general purpose. In this approach, a certain number of interactions are grouped and described by means of a 
multiple scattering theory In condensed simulation (CS), electrons move in finite steps whose length is calculated 
at the beginning of each step, using, for example, the continuum slowing down approximation and a fixed average 
fractional energy loss in the step. 

An alternative approach is provided by the so-called mixed simulation schemes, in which DS is used to simulate the 
£NJ ■ "hard" interactions in which the angular deflections and/or the energy losses are larger than certain cut-off values, 
while the interactions which do not fulfill these conditions, named "soft" , are described within CS approach. 

The statistical models used to determine the energy loss in CS provide good enough results if the tracking steps are 
not too big. On the contrary, the statistical treatment of the angular distribution linked to the multiple scattering 



, is a much more involved task. Since the 40 's a considerable effort has been done to solve the problem and various 
. ^ ■ models (e.g. Goudsmit and Saunderson, Moliere, Lewis) have been developed 0. 

Though multiple scattering theories can describe the angular deflections produced by both elastic and inelastic 
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collisions, it has been a common practice to take care of the effects due to inelastic collisions by correcting the 
distributions calculated considering only elastic collisions Q. Only recently, Negreanu et al. 0] have treated both 
elastic and inelastic collisions on the same footing, using accurate partial-wave differential cross sections. 

The importance of the elastic scattering in the interaction of the electrons with the materials and, in particular, 
k> , in the dosimetry of these particles, has been pointed out by different authors. Andreo et al. [5[ studied in detail the 
limitations of the implementations of the Moliere theory to be used in CS. Li and Rogers Q compared the EGS4 results 
with those obtained by integrating analytically the Moliere distribution, with those quoted in ICRU Report 35 Q 
and with the experimental data of Hanson et al. Q . Urban 0] carried out a comparison of the angular distributions 
obtained with GEANT4 with the same data. 

In the experimental side, the data available are very scarce. To the best of our information, only Kulchitsky and 
Latyshev [lg ] and Hanson et al. [f| have performed measurements of the multiple scattering of electrons in materials. 
Kulchitsky and Latyshev [l(| studied the scattering of 2.25 MeV electrons by foils of different materials from aluminum 
to lead. Their data were in good agreement with the Goudsmit and Saunderson theory for the elements Al, Fe, Cu, 
Mo, Ag and Sn, while for Ta, Au and Pb they found non-negligible differences between theory and experiment. On 
the other hand, Hanson et al. [|[ measured the angular distributions of 15.7 MeV incident electrons scattered by thin 
Be and Au foils. For Au they found discrepancies with the predictions of the Goudsmit and Saunderson theory and 
agreement with the calculations done according to Moliere theory. For Be the experimental data disagreed with the 
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results obtained within the Moliere approach. 

In this work we analyze how different codes describe the scattering of electrons by thin foils of different elements. 
Results for GEANT3, GEANT4, penelope, EGSnrc and MCNPX are compared between them and with the few experi- 
mental data available. In section 2 we give details concerning the simulations we have performed and the MC codes 
which we have used and which are relevant for the calculations performed. Also the way how the data have been 
analyzed is discussed. In section 3 we describe briefly the experiments of Kulchitsky and Latyshev [1(3] and Hanson 
et al. [IJ. Section 4 is devoted to discuss the results we have obtained: the calculations with different versions of 
some of the codes are discussed and the MC results are compared with the experimental data of Hanson et al. Q 
and Kulchitsky and Latyshev [loj and between themselves. In the last section we draw our conclusions. 

II. MONTE CARLO SIMULATIONS 

The simulations have been carried out using an elementary geometry with a monoenergetic pencil beam impinging 
normally on a foil made of a given material and with a given thickness. Both source and foil are in vacuum. The 
incident beam defines the z direction. The directions of the (primary or secondary) electrons emerging from the foils 
define the scattering angles which have been scored in a histogram for angles between and 45 degrees with 180 bins 
of 0.25 degrees. The values are normalized to the solid angle unit given in stereoradian. 

The statistical uncertainty in each bin is calculated as a = »J Q(l — Q)/N with Q the value per history in the 
bin and N the number of histories followed. Through the paper, these uncertainties are given as a number between 
parentheses; e.g., 9.34(1) means 9.34±0.01. 

A. Monte Carlo codes 

As said before, we have performed calculations with various MC codes. All these codes do DS for photons. 
Apart from the possible differences in the photon cross sections, the main differences between the codes are in the 
electron/positron transport. In what follows we quote some details of the codes which are relevant for the MC 
simulations done here. 

1. GEANT3 

GEANT3 [ll[ is a system of detector description and simulation tools designed for high-energy physics at CERN. 
It permits the MC simulation of the transport of elementary particles and different ions in elemental or compound 
materials for energies ranging between 10 keV and 10 TeV. 

The default model for the multiple scattering is based on the Moliere theory, but it is possible to select also the 
pure Gaussian scattering according to the Rossi formulation [l2j]. If the number of Coulomb scatters is below 20, the 
Moliere theory cannot be applied and the simulation is performed within the plural scattering regime in which the 
number of scatters is distributed according to a Poisson distribution. 

The tracking is controlled by means of a series of parameters. DEEMAX is the maximum fraction of kinetic energy 
which a particle can lose in a step. This parameter as well as the multiple scattering introduce upper limits to the 
step length. STEMAX provides an absolute upper limit to the step length, in cm, for each tracking medium. STMIN 
imposes a lower limit to this step length, also in cm, which permits to accelerate the stopping of those particles with 
very small energies. EPSIL determines the boundary crossing precision in cm. In GEANT3 one can use an automatic 
mode for the calculation of these parameter. This is done by selecting AUT0=1. Also, it is possible to fix energy cuts 
for the different particles. In the calculations we have done, these cuts have been fixed to 10 keV for both photons 
and electrons. 

In this work we have used the version 3.2114 of GEANT3. The simulations have been done using the automatic 
mode. To check the feasibility of the corresponding results, we have performed simulations for the Au foil of 18.66 
mg cm -2 with eight different sets of the parameters STEMAX, DEEMAX and STMIN, varying from 0.001 to 0.1, the first 
two, and from 0.0001 to 0.01, the last one. Despite the fact that this is the most exigent case, because correspond to 
the smaller number of interactions, the results obtained arc compatible (within the statistical uncertainty) with those 
provide by the AUT0=1 option. 

2. GEANT4 

GEANT4 code [H} is an object-oriented C++ toolkit which permits the MC simulation of the radiation transport in 
material media for a great variety of particles, materials and energies. For the electromagnetic interactions of photons 
and electrons, GEANT4 permits to use three different physics models: Standard, Low-energy and Penelope. In this 
work the package Low-energy has been used. 
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The multiple scattering approach used in GEANT4 has been developed by Urban [9( based on the Lewis theory [2|]. 
The energy loss is calculated from the actual path length which is computed in every step after performing a path 
length correction. In addition, a lateral displacement and a scattering angle are sampled from given distributions. 

In GEANT4, particles are produced if their energies are above given thresholds which are specified in terms of 
distances, for each volume in the geometry, and internally converted to energy. For the Low-energy model, the 
minimum value for these thresholds is 250 eV. In our calculations we have fixed 10 keV for electrons and 1 keV for 
photons. 

Particles are followed until their kinetic energy is zero, but tracking cuts can be fixed. Other parameters (some 
of them not described in the GEANT4 manuals [Tj, 0j|) provide additional control of electron step. By means of the 
G^UserLimits class it is possible to fix a maximum step size and a maximum track length, as well as the tracking cuts 
above mentioned. The variable dRoverRange determines the maximum fraction of the stopping range that can be 
travelled by an electron in a step. The parameter f inalRange fix the minimum step size below which the electron is 
absorbed locally. Finally, the f r variable permits to control the step size when the electron is transported away from 
a boundary into a new volume. Defaults are dRoverRange=0.2, f inalRange=0.2 mm, fr=0.02 and no maximum 
step size. These are the values we have used in our simulations. 

We have performed test simulations using dRoverRange=0.05, f inalRange=l nm, fr=0.02 and no maximum step 
size. No significant differences have been found with respect to the simulations corresponding to the default values. 

We have used the versions 8.0 (patchOl) of GEANT4 and g4emlow3.0 of the Low-energy package. Some additional 
simulations with the versions 4.1 and 6.0 have also been performed for comparison. 

2. PENELOPE 

penelope [16] is a general purpose MC code which performs simulations of coupled electron-photon transport. It 
can be applied for energies ranging from a few hundred eV up to 1 GeV and for arbitrary materials. Besides, penelope 
permits a good description of the particle transport at the interfaces and presents a more accurate description of the 
electron transport at low energies in comparison to other general purpose MC codes. Details about the physical 
processes considered can be found in [iH ]. 

In penelope electrons and positrons are simulated by means of a mixed scheme where, as said above, collisions 
are classified as "hard" and "soft" . The electron tracking is controlled by means of four parameters. C\ and C2 refer 
to elastic collisions. C\ gives the average angular deflection due to a elastic hard collision and to the soft collisions 
previous to it. C2 represents the maximum value permitted for the average fractional energy loss in a step. On the 
other hand, W cc and W CT are energy cutoffs to distinguish hard and soft events. Thus, the inelastic electron collisions 
with energy loss W < W cc and the emission of bremsstrahlung photons with energy W < W CI are considered in the 
simulation as soft interactions. The maximum step size can be controlled using the parameter s max - 

We have used two version of this code: 2003 and 2005. In what refers to the multiple scattering the main difference 
between both versions is that elastic electron collisions are simulated by means of a model based on the Wentzel 
angular distribution [l6| in the version 2003, while in the 2005 these collisions are simulated by using relativistic 
(Dirac) partial- wave differential cross sections generated by using the computer code elsepa [TtJ ■ 

Some of the simulations performed with penelope have been done with a set of parameters which have been used 
in different simulations in which thin slabs are present (see e.g. [H[). We label these simulations as "safe" and the 
values of the parameters used are: W cc = 5 keV, W cr = 1 keV, C\ = C2 = 0.05. Photons were simulated down 
to 10 keV. Electrons and positrons were absorbed when they slow down to kinetic energies of 100 keV. The mixed 
scheme of the simulation in PENELOPE permits to perform fully "detailed" simulations. We have done calculations 
also in this approach by selecting W cc = 0, W a = —1000 keV, C\ — C2 — 0. In this case the absorption energies were 
fixed to 100 eV for all particles. In all the simulations done with penelope s max was taken to be 1/20 of the width 
of the foil, as it is recommended in [161 ] . 

2. EGSnrc 

EGSnrc [lj| is a general purpose package designed for the Monte Carlo simulation of the coupled transport of 
electrons and photons in arbitrary geometries. Particle energies above a few keV up to several hundreds of GeV can 
be considered. 

Multiple scattering of charged particles is described by means of an approach developed by Kawrakow and Bielajew 
[20l | in which most of the shortcomings of the Moliere multiple scattering theory were fixed. Within this approach, 
track steps are simulated using a single scattering model for short steps and a multiple scattering model for large 
steps. In addition one can select between Rutherford scattering or scattering including both relativistic and spin 
effects. The electron transport algorithm is due to Kawrakow and Bielajew 21] and it is usually known as presta-ii. 

In EGSnrc the simulation is controlled by the following parameters: SMAXIR, which defines upper limits on the step 
size in each region in the geometry; ESTEPR, which fixes the maximum fractional energy loss per electron step in each 
region; ESTEPE, which is a global energy loss constraint, and XIMAX, which gives the maximum first Goudsmit and 
Saunderson moment per step. In our simulations, we have used presta-ii (which is the default transport algorithm) 
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and, as indicated in the manual 0, the default values SMAXIR=10 10 cm, ESTEPR=1, ESTEPE=0.25 and XIMAX=0.5 
must not be changed. The version V4-r2-2-3 has been used for calculations. Photons and electrons were followed 
down to 10 keV. 

2. MCNPX 

MCNPX [22j is a general purpose MC code which permits the description of the transport of different particles in 
arbitrary materials. Photons, electrons and neutrons, as well as other 29 particles between leptons, baryons, mesons 
and even light ions can be considered. The upper energy limits for electrons and photons are 1 and 100 GeV, 
respectively. A lower limit of 1 keV is fixed for these particles. 

The angular deflections in the multiple scattering of electrons are calculated according to the Goudsmit and Saun- 
derson theory. When electron energies are below 0.256 MeV, the corresponding cross sections are obtained from 
numerical tabulations developed by Riley et al. (23j . For higher energy electrons, the cross sections are approximated 
as a combination of the Mott and Rutherford cross sections, including a correction factor which takes care of the 
screening. 

The particle tracking is governed by EMCPF, which is the upper energy limit for detailed photon physics treatment, 
EMAX, which fixes the upper limit for electron energy, and the low energy cutoffs. We have used the default value (100 
MeV) for EMCPF and EMAX=20 MeV, while the low energy cutoffs were fixed to 1 keV. In our simulations we have used 
the version 2.4.0. 

In the calculations done with MCNPX, the angular distributions obtained by using the binning indicated in Sect. [TT1 
showed a stepwise shape in which various neighbour bins had the same value. This is due to the fact that the Goudsmit 
and Saunderson angular distribution is sampled by using the inverse transform method applied to a histogram of 34 
bins with widths equal or larger than 1 degree. In the cases analyzed here, where few interactions occur, this lack of 
precision in the sampling produces wrong results. To solve this problem, the bins in the simulations performed with 
MCNPX have been enlarged to 1 degree. 

III. DESCRIPTION OF THE EXPERIMENTS 

In the two experiments mentioned above, the angular distribution of electrons scattered by foils of different materials 
and widths were measured. 

Kulchitsky and Latyshev [l(| studied the scattering of 2.25 MeV electrons by foils of Al, Fe, Cu, Mo, Ag, Sn, Ta, Au 
and Pb. The thicknesses (see first row in Table HJ were chosen in order to maintain the half- width of the distribution 
~10 degrees in all cases, thus ensuring the validity of the multiple scattering theory. Measurements for angles between 
and 35-40 degrees were done. Electrons with energies in the range 2250±23 keV formed their beam. 

The data obtained were compared to the Goudsmit and Saunderson multiple scattering theory predictions. Specifi- 
cally, the comparison was done for the angle at which the angular distribution reduces by a factor 1/e. Kulchitsky 
and Latyshev estimated the uncertainty in 0\/ e as 3-4% at most, but they provided only relative values of the distri- 
butions measured, what impedes absolute comparisons with our calculations. 

Hanson et al. Q measured the angular distributions of electrons scattered by thin Be (257 and 495 mg cm -2 ) and 
Au (18.66 and 37.28 mg cm -2 ) foils for 15.7 MeV incident electrons. In the case of Be foils, the measurements were 
done up to 6 degrees. For Au foils the angular range was extended up to 30 degrees. The beams included electrons 
with energies within 6 percent of the maximum energy. 

Hanson et al. compared their experimental data to the predictions of the Goudsmit and Saunderson and the Moliere 
theories for the angle 9\/ e . They do not quote the uncertainties of their measurements but provide absolute values of 
the distributions, what allows a more complete comparison. 

The results obtained with the different MC codes have been compared to the experimental data by means of the 
quantities Q\j e and, when possible, F(0). This last corresponds to the maximum of the angular distribution, at 6 = 0. 





Al 


Fe 


Cu 


Mo 


Ag 


Ta 


Au 


Pb 


thickness [mg cm 2 ] 


26.6 


15.4 


17.15 


12.4 


11.55 


8.9 


8.9 


7.9 


Of/1 [deg] 
&l7e t de g] 


9.50 
9.49 


9.60 
9.33 


10.40 
10.61 


10.25 
9.97 


10.20 
10.21 


9.85 
9.91 


9.20 
10.07 


9.70 
9.36 



TABLE I: Comparison of the values of the parameter obtained after reanalyzing the experimental data of Kulchitsky and 
Latyshev [l(j with those quoted by these authors. The thicknesses of the different foils are also given. 
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0i H /T [deg] 


KTe t de s 


F Han (0) [sr- 1 ] 


F cxp (0) [sr- 1 ] 


Be (257 mg 


cm- 2 ) 


3.06 




104.06 




Be (495 mg 


cm- 2 ) 


4.25 




50.23 




Au (18.66 m; 


; cm -2 ) 


2.58 


2.50 


144.77 


143.29 


Au (37.28 mj 


5 cm -2 ) 


3.76 


3.71 


65.66 


64.68 



TABLE II: Comparison of the values of the parameters #™ p and F exp (0) obtained after reanalyzing the experimental data of 
Hanson et al. [§| with those quoted by these authors. We have not calculated the values for the Be foils because the angular 
distributions are not available in Ref. 



4.1 6.0 8.0 


Be 


F(0) 


[sr" 1 ] 


66.13(8) 40.63(4) 41.84(5) 


(495 mg cm -2 ) 




[deg] 


3.762(4) 4.778(5) 4.735(5) 


Al 


F(0) 


[sr- 1 ] 


16.23(2) 11.44(1) 10.83(1) 


(26.6 mg cm - ] 


6l/e 


[deg] 


7.473(9) 9.17(1) 9.33(1) 


Ta 


F(0) 


[sr" 1 ] 


10.13(1) 9.72(1) 9.34(1) 


(8.9 mg cm" 2 ) 




[deg] 


9.51(1) 10.03(1) 10.04(1) 



TABLE III: Values of the parameters F(0) and 9\i e obtained for the calculations shown in Fig. [T]and performed with the three 
versions of GEANT4 considered in this work. The uncertainties are given at 1<t level. 



To do that we have fitted a Gaussian function 



F(0) = F(0) exp 



(1) 



to the simulated angular distributions, in the angular range between and the first angle larger than 9\/ e . The 
Levenberg-Marquadt method [24| has been used. 

The experimental data of Kulchitsky and Latyshev and Hanson et al. § have been reanalyzed following the 
same criteria used for the simulated distributions. In this way we have recalculated the "experimental" #w e and F(0). 
We call the values obtained in this way 6^? and F cxp (0), respectively. We give them in Tables HI and [Til where they 

are compared to the values quoted in Refs. [l(|, and fS], an< ^ ^ Han (0). The values F oxp (0) have not been 

calculated for the experiment of Kulchitsky and Latyshev because the absolute distributions are not available in [101 ] . 
On the other hand, the values and F exp (0) have not been calculated for the Be foils in the experiment of Hanson 
et al. because they do not provide the corresponding angular distributions. 

In the case of the experiment of Kulchitsky and Latyshev (see Table [I]), the maximum difference between our results 
and those quoted in Ref. [lfj is for Au foil (8.6%). In the case of the Au foils in the experiment of Hanson et al., the 
maximum difference is smaller than 3.5% (see Table HT|) . 

To finish we mention here the fact that there is a mistake concerning the Sn data in the experiment of Kulchitsky 
and Latyshev. The thickness quoted in Table 1 of Ref. [T3| is incompatible with the value 6\/ e =10.9 degrees they 
obtained. All the simulations we have performed disagree by more than 20% with this value. This disagreement can 
be seen also in Fig. 5 of Ref. [H where the data were compared to calculations using the Lewis multiple scattering 
theory. We have not considered this foil in our discussion. 



IV. RESULTS 



A. Comparison of different code versions 



In the case of GEANT4 and penelope, different versions have been considered in the calculations. Here we compare 
the results obtained between them. 

In Fig. [T]we show the angular distributions obtained with the versions 4.1, 6.0 and 8.0 of GEANT4 for three different 
cases: 15.7 MeV electrons impinging on a Be foil with 257 mg cm -2 (upper panel) and 2.25 MeV electrons incident on 
Al (medium panel) and Ta (lower panel) foils with 26.6 and 8.9 mg cm -2 , respectively. As we can see, the calculations 
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GEANT4 



Be (495 mg cm 2 ) 
B e =15.7 MeV 

v. 8.1 
v. 6.0 
v. 4.1 
i 



20 
15 

ST 

j 


12 

V 

6 

3 



5 10 15 20 25 

[deg] 

FIG. 1: Angular distributions obtained with the versions 8.0 (solid curves), 6.0 (dashed curves) and 4.1 (dotted curves) of the 
GEANT4 code. We show the results for 15.7 MeV electrons impinging in a Be foil with 495 mg cm -2 (upper panel) and for 
2.25 MeV electrons incident on an Al (medium panel) and Ta (lower panel) foils with 26.6 and with 8.9 mg cm -2 , respectively. 
Statistical uncertainties (la) are smaller than 2.5% and have not been plotted. 

performed with the versions 6.0 and 8.0 provide rather similar results, while those found for the 4.1 version disagree 
with them. The differences are larger the smaller the atomic number of the foil material is. In order to quantify 
these differences we show, in Table UTT1 the values of the parameters F(0) and 9i/ e corresponding to these angular 
distributions. For Be and Al, 8i/ e is 25% larger for the version 8.0 than for the 4.1 one, while for Ta is less than 6%. 
On the other hand, for Be and Al, the F(0) value obtained with the version 4.1 is more than 50% larger than for the 
version 8.0, while is only 8% larger for Ta. From these results it is evident that the various changes in the details 
concerning the multiple scattering theory carried out from version 4.1 till version 8.0 (see Q) are not at all negligible. 

For penelope, the versions 2003 and 2005 have been considered and for each one, both "detailed" and "safe" 
simulations have been performed. In Fig. [2] we show the corresponding angular distributions. Full and dotted lines 
represent the results of detailed calculations with versions 2005 and 2003 respectively. Dashed and dashed-dotted 
curves correspond to the "safe" simulations. The values of the characteristic parameters are given in Table ITVl Two 
findings deserve a comment. First, the results obtained with the "safe" simulations are very similar to those obtained 
with the "detailed" simulation for the two Be foils (left panels). In these cases the differences are 2% at most. 
However, for the Au slabs, the differences between the results obtained in both simulations are larger. For F(0), 
these differences are above 20% for the version 2003 and around 15% for the 2005. The differences between the values 
obtained for #i/ e within "safe" and detailed simulations are around 5% for the version 2003 and slightly larger for 
the version 2005. The discrepancies observed between "safe" and "detailed" simulations are linked to the number of 
interactions suffered by the electrons in the foils. In case of Be, electrons suffer ~ 10 hard interactions in average; for 
Au foils, the average number of hard interactions is smaller than 1. As indicated in [lj|, this makes "safe" simulations 
to be into agreement with the "detailed" ones for Be foils and to show the differences discussed in the case of the Au 
foils. By fixing adequately the tracking parameters these differences can be strongly reduced. For example, in the 
case of the thin Au foil, a reduction of the parameter C\ to 0.01 makes the difference at 9 = degrees to diminish 
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PENELOPE 

160 



120 



40 





v. 2005 


~.\\ 


detailed 








O Hanson et at 




\ Au (18.66 mg cm" 2 ) 




1 1 ^>-L^J^_ | „ 




FIG. 2: Angular distributions obtained with penelope for the four foils of the experiment of Hanson et al. |8(. Detailed 
simulations with the versions 2003 (dotted curves) and 2005 (solid curves) and simulations performed with the "safe" parameter 
indicated in the text also with versions 2003 (dashed-dotted curves) and 2005 (dashed curves) are shown. For the two Au foils, 
the open circles represent the experimental data. Statistical uncertainties (la) are smaller than 2% and have not been plotted. 









2003 2005 








"safe" detailed "safe" detailed 


Be 

(257 mg cm" 2 ) 


F(0) 


[ST' 1 ] 

[deg] 


75.9(1) 76.5(1) 95.8(2) 97.0(2) 
3.640(6) 3.620(6) 3.173(5) 3.141(5) 


Be 

(495 mg cm" 2 ) 


F(0) 


[sr- 1 ] 
[deg] 


37.51(8) 37.73(8) 45.7(1) 46.1(1) 
5.18(1) 5.15(1) 4.595(9) 4.563(9) 


Au 

(18.66 mg cm -2 ) 


F(0) 


[sr" 1 ] 
[deg] 


107.8(2) 119.3(3) 126.6(3) 146.9(3) 
3.033(5) 2.796(5) 2.791(5) 2.480(4) 


Au 

(37.28 mg cm -2 ) 


F(Q) 


[sr" 1 ] 
[deg] 


51.0(1) 53.9(1) 59.4(1) 64.4(1) 
4.376(8) 4.181(8) 4.032(8) 3.763(7) 



TABLE IV: Values of the parameters F(0) and 6i/ e obtained for the calculations shown in Fig. [2] and performed with the two 
versions of penelope considered in this work. The uncertainties are given at la level. 



from 14.7 to 6.7%. 

The second point to be noted concerns the differences observed between the two versions of penelope. By com- 
paring the results corresponding to the detailed simulations, differences between 19 and 27% (for F(0)) and 11 and 
15% (for #i/ e ) are found, the larger differences occurring for the thin Be slab and the smaller for the thick Au foil. 

In addition, the substitution of the modified Wentzcl angular distribution (used in the version 2003) by the relativis- 
tic partial-wave differential cross sections (which incorporates the version 2005) permits, in the case of the detailed 
simulations, a very good description of the experimental data for the Au foils (which are plotted in Fig. [5] with open 
circles). 

In order to go deeper into the differences observed between the two versions of penelope, we have calculated 
the mean values (/z) and (/j?) of the quantity /i = (1 — cos#)/2 which is commonly used to measure polar angular 
deflections produced by single scattering events, instead of the scattering angle 9. The results of the calculations 
for the material foils considered in the experiment of Kulchitsky and Latyshev [Tc| are shown in Fig. [3J In these 
calculations both the forward (0 < 6 < 90 deg) and the backward (90 deg < 8 < 180 deg) distributions have been 
considered to perform the corresponding integrals. It should be pointed out the almost perfect agreement between 
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the results obtained with the two versions of the code. This was expected according to the fact that, as indicated in 
the PENELOPE (v. 2003) manual [l6], the multiple scattering theory used in this version (the modified Wentzel model) 
was fixed to reproduce the values of (/z) and (/i 2 ) obtained for the actual partial- wave differential cross section which 
are used, instead, in Penelope (v. 2005). In principle, this equality ensures a reasonable description of the multiple 
scattering in "normal" simulation conditions with enough interactions. It is clear that in the cases analyzed here 
(where electrons go through rather thin foils) the differences between the models used show up. 

In what follows we have considered only the version 8.0 of GEANT4 and the detailed simulations of the versions 
2003 and 2005 of penelope. 

B. Comparison with experiment 

First we have compared the results of the simulations to the experimental data of Kulchitsky and Latyshev 
In Fig. [4] we show the ratio of the values 6\/ e , obtained with the various MC codes, to the values 0y^, obtained in 
our reanalysis of the corresponding experimental data, as a function of the atomic number Z of the foil materials. In 
addition, the ratio Oy^/ff^? is plotted as reference (see circled plus symbols.) 

While EGSnrc, GEANT3 and penelope (v.2005) underestimate in all cases the values of 0f^, penelope (v.2003) 
overestimates it for all foils and MCNPX and GEANT4 show, in this respect, a Z dependent behaviour. The differences 
between the values Q\u obtained for the different codes and those we have found in our reanalysis of the experimental 
data, which show the maximum and minimum absolute values, are shown in Table [V] As we can see, EGSnrc and 
PENELOPE (v. 2005) maintain these differences within a range of 3%, GEANT3, GEANT4 and PENELOPE (v. 2003), within 
8%, while mcnpx shows a range bigger than 14%. 

In order to quantify globally these differences we have calculated, for each MC code, the sum 



/e u l/c 



eiTe) , (2) 




E e =2.25 MeV 

















FIG. 3: Values of (fx) (upper panel) and (/i 2 ) (lower panel) as a function of the atomic number Z of the foil materials in the 
experiment of Kulchitsky and Latyshev [10] . Detailed simulations with the versions 2003 (open circles) and 2005 (black points) 
of penelope are shown. Statistical uncertainties (la) are smaller than the symbols used to plot the results. 
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□ EGSnrc © Kulchitsky & Latyshev 
o MCNPX 
A GEANT3 




0.9 - 

A 
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E =2.25 MeV 



20 40 60 80 100 

Z 

FIG. 4: Ratio of the values obtained with the various MC codes to the values 9^^, obtained in our reanalysis of the 
corresponding experimental data as a function of the atomic number Z of the foil materials in the experiment of Kulchitsky 
and Latyshev [lQ|| . Circled plus symbols correspond to the ratio between the values quoted by Kulchitsky and Latyshev 
(see Table U]) and SfJJ. The lines between symbols are only to guide eyes. Uncertainties are considered at the la level and in 
most cases are smaller than the symbols used to represent the data. 



EGSnrc 


MCNPX 


GEANT3 


GEANT4 


pen 2003 


pen 2005 


-8.7% (Ag) 
-6.5% (Fe) 


10.9% (Mo) 
-3.7% (Cu) 


-11.2% (Ta) 
-3.5% (Al) 


4.4% (Pb) 

-3.5% (Cu) 


11.6% (Pb) 
4% (Ag) 


-7.8% (Cu) 
-5.1% (Pb) 


S [deg 2 ] 4.9 


1.9 


5.2 


0.5 


4.3 


3.1 



TABLE V: Relative differences (#i/ e — #i/e)/$i/e with maximum (first row) and minimum (second row) absolute values, obtained 
for each MC code for the foils in the experiment of Kulchitsky and Latyshev [ljj. The corresponding foil materials are given. 
Last row shows the values of S as defined by Eq. @. 

where i runs over the eight different foils considered. The better codes correspond to those with the smaller values of 
S. The results obtained are shown in the last row of Table fVl On the base of this index, the "best" results are those 
obtained with GEANT4 while the large discrepancies are obtained with GEANT3, EGSnrc and Penelope (v. 2003). 

A slightly different situation is found when we compare the results of our simulations to the experimental data of 
Hanson et al. 8]. In Fig. Owe show the ratio of the values 0i/ e (upper panels) and F(0) (lower panels) obtained with 
the various MC codes to the values and F cxp (0), respectively, for the two Au foils in the experiment of Hanson 

et al. 0. The dashed lines represent the ratios between the values O^ff and F Han (0) quoted by Hanson et al. (see 
Table |TT| and or F c * p . As we can see, EGSnrc, GEANT3, GEANT4 and penelope (v.2005) provide a very good 
description of the experiment, mcnpx slightly overestimates (underestimates) O^P (F exp (0)) and penelope (v. 2003) 
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Au (18.66 mg cm 2 ) 



Au (37.28 mg cm" 




I T A---' 



E. =15.7 MeV 



FIG. 5: Ratio of the values 0\/ e (upper panels) and F(0) (lower panels) obtained with the various MC codes to the values 
and 7 ?exp (0), respectively, obtained in our reanalysis of the data for the two Au foils in the experiment of Hanson et al. 
The dashed lines represent the ratios between the values 9yf and i ?Han (0) quoted by Hanson et al. (see Table Hi]) and 
F cxp . Uncertainties are considered at the lcr level and are smaller than the symbols used to represent the data. 
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FIG. 6: Ratio of the values 0\/ e (upper panels) and F(0) (lower panels) obtained with the various MC codes to the values 
0^J£ and J ?Han (0), respectively, quoted for the Be foils in the experiment of Hanson et al. || (see Table HI)) . Uncertainties are 
considered at the la level and are smaller than the symbols used to represent the data. 



gives the "worst" results. The excellent agreement found for GEANT4 is not surprising because the multiple scattering 
theory developed by Urban and which is used in this code has been tuned to reproduce these Au experimental data 
1- 

Fig. [5] shows the comparison with the experimental results of Hanson et al. [8j for the two Be foils. As said before, 
data are not available for these foils and we could not perform the corresponding reanalysis. In this case we have 
used directly the values and F Han (0) quoted by Hanson et al. (see Table |H]) as reference values. Here the "best" 
description of the experiment is provided by EGSnrc, MCNPX and PENELOPE (v. 2005). GEANT3 and GEANT4 show very 
similar results and the major differences occur for penelope (v. 2003), reaching values around 20% or even larger, 
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both for 9\u and F(0). On the other hand, and as in the case of the Au slabs, relative results are rather similar for 
both foils. 
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FIG. 7: Angular distributions obtained with the various MC codes for the four foils of the experiment of Hanson el al. 
Statistical uncertainties (la) are smaller than 2% and have not been plotted. 



C. Comparison of the simulations 



To finish the discussion we compare now the angular distributions obtained with the different simulation codes. 

In Fig. [7] we show these distributions for the four foils of the experiment of Hanson et al. [8]. As we can see, the 
behaviour observed is different in the case of the Be foils (left panels) than in the Au ones (right panels). For Be, 
the calculations performed with EGSnrc (dashed-dotted curves) show the maximum values F (0), and are in rather 
good agreement with the results obtained with Penelope (v. 2005) (solid curves) and MCNPX (solid triangles). The 
GEANT3 (dashed curves) and GEANT4 (open squares) results overlap for all values of 9 and, finally, the PENELOPE 
(v. 2003) (dotted curves) are clearly below the other calculations at small angles. The relative differences between the 
maximum and minimum values of F(0) are 33% and 26.3% for the thin and the ticker foils, respectively. 

In the case of the Au foils (see right panels), the situation of the Penelope (v. 2003) results is the same, but the 
other calculations are all of them grouped. This is particularly clear in the case of the thicker foil (lower right panel). 

The results obtained for the different foils considered by Kulchitsky and Latyshev [l(| show up a different behaviour. 
In Fig. [8] we have plotted the angular distributions corresponding to the Al, Fe, Ag and Au foils. Similar results 
are found for the other four foils studied in this experiment, penelope (v. 2003) provides the smaller F(0) values in 
all cases except for Mo, in which MCNPX is the code showing the minimum. EGSnrc shows the larger values for Al, 
Fe and Cu and GEANT3 for the remaining foils. As we can see, the results appear to be spread for forward angles 
and the spreading increases with the atomic number of the material. The relative differences between the maximum 
and minimum values of F(0) range between 17.3% for Cu to 45.9% for Pb, being ~ 20% for Al and Fe. In all cases 
analyzed, penelope (v. 2005) and EGSnrc are in good agreement. 

To finish, it is worth to point out the situation observed for the Au foils. As we can see (lower right panel in Fig. [8]) 
the various codes produce rather different angular distributions for 9 < 8 degrees. This contrasts with the overlapping 
shown by the results corresponding to the Au foils in Fig. [7J (see right panels.) 
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FIG. 8: Same as in Fig. [7] but for the Al, Fe, Ag and Au foils in the experiment of Kulchitsky and Latyshev [ic| . 

V. CONCLUSIONS 

In this work the experimental data available for multiple scattering by thin material foils have compared to the 
simulation results obtained with the MC codes EGSnrc, mcnpx, GEANT3, GEANT4, penelope (v.2003) and PENELOPE 
(v.2005). 

Simulations performed with various older versions of GEANT4 show up non-negligible differences in the angular 
distributions with respect to the version 8.0, differences which are bigger the lower the Z of the material is. 

Different simulations done with penelope have indicated large differences in the angular distributions between 
the versions 2003 and 2005, though the mean values (//) and (/j, 2 ) obtained with both versions are very similar, as 
expected. In addition, we have found discrepancies between the "safe" and "detailed" simulations which are related 
to the number of interactions suffered by the electrons in the foil and which can be largely reduced by an adequate 
selection of the tracking parameters. 

We have compared the results of our simulations with the experimental data of Kulchitsky and Latyshev [lol ] by 
means of the characteristic angle 1 / e . We have found that both EGSnrc and penelope (v. 2005) show differences 
with respect to c *^ e which are within a range of 3% for all foils. However, the better description of the experimental 
data is provided, globally, by GEANT4. 

The comparison with the results of Hanson et al. [1] shows that EGSnrc, GEANT3, GEANT4 and penelope (v.2005) 
give a good description of the data in the case of the Au foils. For the Be foils, the better agreement with the 
experiment occurs for EGSnrc, mcnpx and penelope (v.2005). 

All codes (except penelope (v. 2003)) show a good agreement for the two Au foils in the experiment of Hanson et 
al. For the Be foils, the maximum relative differences between the codes at 9 — are larger than 25%. In the case 
of the foils in the experiment of Kulchitsky and Latyshev, these maximum relative differences grow with Z , reaching 
~ 50% in the case of the Pb foil. In general, EGSnrc and penelope (v.2005) are in rather good agreement for all 
cases studied. 

The effects observed in these thin material foils appear to be non-negligible. Differences between the various 
MC codes are relevant and the propagation of them to other situations should be investigated in detail. One of 
these situations concerns the simulation of electron linear accelerators, where different thin foils are present. The 
combined effect of all of them could modify the conclusions found in previous works in which various codes have been 
intercompared (e.g. in Ref. fl8l|.) 

On the other hand it should be interesting to perform new experiments with different materials in order to permit 
a complete test of the multiple scattering theories used in the MC simulation codes of the radiation transport. 
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